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A five dimensional brane cosmology with non-minimally coupled scalar field to gravity is con- 
sidered in the Jordan frame. We derive the effective four dimensional field equations on a 3 + 1 
dimensional brane where the fifth dimension is assumed to have an orbifold symmetry. We show 
that energy conservation still holds for matter on the brane in spite of the existence of a scalar 
field (the Brans-Dicke field) in the bulk. Finally, we discuss some cosmological consequences of this 
setup. 

o: 

O , I. INTRODUCTION 

O; 

^-^ i Recently the possibility that the observable Universe is a brane-world m embedded in a higher dimensional space- 

t> time has generated a great deal of interest. Such a claim is motivated by the strongly coupled sector of E^ x E^ 

O , heterotic string theory which can be described by a field theory living in 11 dimensional space-time |0]. The 11 

dimensional world is comprised of two 10 dimensional hypersurfaces embedded on the fixed points of an orbifold and 

the matter fields are assumed to be confined to these hypersurfaces which in this scenario are known to be 9 branes. 

fvj ' After compactifying the 11 dimensional theory on a Calabi-Yau three fold, one obtains an effective 5 dimensional 

theory which has a structure of two 3 branes situated on the orbifold boundaries. The theory allows A^ = 1 

C^ ' supergravity with gauge and chiral multiplets on the two 3 branes. Motivated by the above setup there has been an 

^ ' attempt to understand the case where the bulk is 5 dimensional anti de-Sitter space. It was shown that gravity could 

t^ , also be localized in this scenario Q . These authors demonstrated that in a background of a special non-factorizable 

geometry an exponential warp factor multiplies the Poincare invariant 3-1-1 dimensions in the metric. The model 

consists of two 3 branes situated rigidly along the 5th dimension compactified on a S^ /Z2 orbifold symmetry. In order 

^~^ ' to realize this simple scenario the two branes must have opposite tensions. 

f^ \ It was realized that the brane-world construction could modify the early Universe considerably and, furthermore, 

^^ i the two branes with opposite tensions should behave differently ||5[] . The most important observation was the departure 

^j ' of the evolution equation on the brane from the standard four dimensional evolution equation when no branes are 

Q^, present |g] . The presence of branes and the requirement that the fields should be localized lead to a non-conventional 

J ■ cosmology which requires a more exhaustive study. Some attention has been devoted to the effective gravity induced 

(^ on the brane H and a great emphasis was placed on inflationary cosmology B and, more recently, post-inflationary 

IL" ■ brane cosmology has been considered in Ref. |^ . Apart from inflationary and post-inflationary brane cosmology, there 

. ^ \ has been a great deal of interest in the four dimensional cosmological constant problem |lfl-0] . These authors have 

j^ • considered a five dimensional action with a scalar field non-minimally coupled to the five dimensional gravity and to 

^ ' the four dimensional brane tension. There has also been some discussion on the localization of gravity in this setup 

|l3| . In this paper we would like to extend the calculation made in Ref. [|| to include the dynamics of a scalar field 

which lives on the brane and the bulk. We then study the effective field equations on the 3 -f 1 dimensional brane 

locally and globally in presence of the other brane which is assumed to be rigidly located on the orbifold symmetry 

along the 5th dimension. Finally, we discuss the pertaining cosmology. 

II. BRANS-DICKE BRANE MODEL IN 5 DIMENSIONS 

In this article we restrict ourselves to one extra space-like dimension with a scalar field non-minimally coupled 
to gravity in 5 dimensions. We analyze the physics in the Jordan frame where the interpretation of the results is 
straightforward, as we shall see in the next sections. We believe that such a scalar field in the five dimensional theory 
can be a dilaton which is purely an outcome of dimensional reduction from some higher dimensional theory to 5 
dimensional space-time n. We assume the following action for the rest of our discussion 
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where R is the Ricci scalar associated with the S-dimensional space-time metric ^ab , <^ is a scalar field which we will 
call Brans-Dicke (BD) field, w is a dimensionless coupling constant which determines the coupling between gravity and 
the BD scalar field and finally, £„! represents the Lagrangian for the matter fields. Latin indices denote 5 dimensional 
components (A,B = 0, . . . , 5) and for our convenience we choose k| = 1. The variational derivative of the action 
Eq. (m with respect to ^ab and (j) yields the field equations M^ 
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where T — T q is the trace of the energy-momentum tensor of the matter content of 5 dimensional space-time. 
Notice the Zoj + A denominator in the right hand side of the BD field equation instead of the familiar 2a; -I- 3 in the 
4-dimensional casejl^]. This is determined by requiring the validity of the equivalence principle in our setup (see 
Weinberg's book |15] for a discussion of this topic in the context of 4-dimensional Brans-Dicke theory). 

Before we discuss the energy-momentum tensor, let us define the five dimensional metric which has the following 
form 
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where ^,1^ ~ 0, . . . , 3, and y is the coordinate associated with the fifth dimension which is assumed to be compact 
with a range —1/2 < y < 1/2. We also assume an orbifold symmetry along the fifth direction y ~^ ~y- As we shall 
see in the coming sections, this will help us to simplify our calculation. Since the opposite points along the fifth 
dimension are identified, we will only be interested in the interval < y < 1/2, which is the spacing between the two 
branes situated a,i y — Q and y — 1/2. Next we define the energy-momentum tensor 



T^B — T Blbulk + T" Bibranc: 



(5) 



where the subscripts brane and bulk suggest the corresponding energy-momentum tensors. For simplicity we assume 
the bulk is devoid of matter other than the BD scalar field. The brane matter fields are stuck at y = and y = 1/2 
with the following energy momentum tensor 
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where the subscript '*' corresponds to the matter on the brane at y = 1/2. While writing the above expressions we 
have taken the infinitely thin brane limit. 

III. EQUATIONS OF MOTION 

Since we are interested in exploring the flat cosmology, we consider a 5 dimensional flat metric ansatz of the following 
form 



ds^ — —n^{T, y)dT'^ + cl^{t, y)Sijdx^dx^ + b^i^r, y)dy^ , 



(8) 



where i,j = 1,2,3. With this metric ansatz we are now able to write the equations of motion. The (0,0) component 
reads 
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the {i,j) component is given by 
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the (0, 5) component takes the form 



and, finally, for the (5, 5) component one has 
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The equation of motion for the BD field reads 
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where the dot corresponds to the time derivative with respect to r and the prime corresponds to derivatives with 
respect to y. We make the assmxiption that the metric and the BD field are continuous across the branes localized 
at y = and y — 1/2. However, their derivatives can be discontinuous at the brane positions in the y direction. 
This suggests the second derivatives of the scale factor and the BD field will have a Dirac delta function associated 
with the positions of the branes. Since the matter is localized on the branes it will introduce a delta function in the 
Einstein equations which will be matched by the distributional part of the second derivatives of the scale factor and 
BD field. For instance at j/ = 0, we have 
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where the hat marks the non-distributional part of the second derivative of the quantity. The part associated with a 
delta function, [/'], is a jump in the derivative of /. Here / could be any of the three quantities a, n or (jj. The jump 
in / at y = can be written as 



[/']-/'(o+)-/'(o-) 

and the mean value of the function / at y = is defined by 
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After substituting Eq. (14) in the Einstein equations it is possible to find out the jump conditions for a and n by 
matching the Dirac delta functions appearing in the left-hand side of the Einstein equations to the ones coming from 
the energy- momentum tensor Eqs. (P^R)- For the BD field one has to use Eq. (O) to evaluate the jump condition. 
We get the following results for the jump conditions 
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where the subscript '0' stands for the brane at y = 0. Similar conditions hold for the brane fixed at y = 1/2 by 
replacing '0' by '1/2' and p and p by p* and p,. The first two conditions, Eqs. (|l^) and (|l8|), are equivalent to 
Israel's junction conditions in general relativity ||l^ (see Ref. [|| for a discussion of its application in the context of 
brane world). It is important to note that the above jump conditions aX y — Q depend on the energy density and the 
pressure component of the brane world. Remarkably for the radiation dominated phase on the brane, p — 3p, the 
jump condition for (f) vanishes suggesting that the BD field takes a constant value on the bulk and the brane. This 
has important consequences which we shall discuss later on. 

Taking the jump of the (0, 5) component of the Einstein equation and substituting Eqs. ( [l7Hl8| ) we get the continuity 
equation for the matter on the brane 
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The fact that the energy content of the brane is still conserved in this scenario seems to be at odds with recent results 
obtained independently by several authors who conclude that the presence of a dilaton in the bulk will lead to a non- 
trivial coupling with the matter on the brane which from the point of view of an observer living on the brane would 
be seen as matter leaking from the brane p^j2^ -|2^ . The reason why our situation is different is because the coupling 
between the BD field and ordinary matter on the brane was chosen in order to satisfy the equivalence principle, as 
was pointed out above. Had we chosen the coupling in (n3)to be the usual 4-dimensional value (2a; + 3)~^ we would 
end up with a situation where the conservation equation^ 29) would not hold and energy could leak from the brane. 
Taking the mean value (in the sense of Eq. (ttq)) of the (5, 5) component of Einstein's equations we can now obtain 
a Friedmann type equation on the brane by following a very similar procedure to the one introduced in [p|. Using 
the fact that due to the orbifold symmetry y <-> — y, jj/'ft = (for / any of the quantities a, n or 0), we can discard 
all the terms involving mean values in the average of the (5, 5) component of the Einstein equations. The equation so 
obtained will still involve a term containing 0, but we can use the mean value of the BD field equation (13) to write 



this in terms of a, h and n, and their derivatives. After a tedious calculation we obtain the Friedmann type equation 
on the brane 
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While deriving the above equation we have assumed the extra dimension to be static, b = bo, from the point of view 
of the brane observer and we have also fixed time in such a way that no = 1. This corresponds to the usual choice 
of time in conventional cosmology. It is interesting to note that by taking the limit w — *■ cxd in the right-hand side of 
Eq. (pl|), we obtain exactly the same expression as in general relativity (compare with Eq. (20) in Ref. Q). However, 
the above expression differs from the corresponding equation of motion without the presence of branes (see Ref. |Q6[ ) , 
namely due to the presence of quadratic terms in p and p and the absence of linear terms in the energy density and 
pressure in ([2l|). 

From the mean value of the BD field equation we obtain an equation of motion for (j) on the brane 
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Note that in order to obtain Eq. ( p2D we also have to assume the non-distributional part of (p" vanishes, otherwise, 
a term involving 0" will appear in the BD field equation. As we shall see in the next Section it is possible to 
obtain cosmologically interesting solutions which verify this condition. Together with the conservation equation (pO), 
Eqs. ( pi| ) and ( p^ determine the cosmology of the brane at y = 0. Comparing equation ( |2^ ) with the usual BD field 
equation in 4 dimensional BD theory with no branes we note there is an extra power of p — 3p and an extra factor of 
0Q ^ on the right hand side of the BD field equation in the brane model. 



IV. TOPOLOGICAL CONSTRAINTS 



As in the general relativity case, in our setup matter in one brane is also constrained by matter on the second brane 
(see Ref. Q for a detailed discussion of this point). In order to get the constraints, we proceed as in Ref. [^. We will 
use the following ansatz for the solution 
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where the subscripts '0' and '1/2' denote the branes at y = and y = 1/2 with a similar expression for n. As to the 
BD field 0, we assume it has the following form 

^-Mt)+^^\y\- (24) 

Notice that the term quadratic in y which is present in a and n is absent from (f>. This difference will be analyzed 
below. Since no(t) is a completely arbitrary function which can be fixed by a choice of time coordinate we fix no(i) = 1 
in order to get the conventional definition of time. For b we use 

b^bo , (25) 

where bo is a constant. This ansatz corresponds to stabilizing the fifth dimension. This is by no means a trivial 
assumption but for simplicity we assume the fifth dimension to be static. Using the above ansatz for the static extra 
dimension the (0, 0) component of the Einstein equations for the brane situated at y = yields 
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In order to be consistent with Eq. ( pi[ ) the right-hand side of this equation must be quadratic in the energy density. 
Noticing that the jump conditions for a and n are linear in the energy density, this means that the right-hand side of 
Eq. (Um must contain only terms quadratic in the jump for the various quantities involved. Terms linear in the jump 
must therefore vanish yielding the condition 

[a']o = -[a'] 1/2 . (27) 

Following a similar procedure with the {i,j) components of the Einstein equations, with the help of the previous 
constraint Eq. ( p7| ) we obtain our last constraint 

Mo = -Ml/2. (28) 

These constraints have the same form in terms of the jumps in a and n as in the general relativity case. However, 
using the jump conditions Eqs. ( |l7|) and ( |l8|) they result in different relations between the energy density and pressure 
in both branes than the ones found in Ref. [pj. 

Assuming the matter on the branes obeys an equation of state p ~ ep for the brane at y = and p^ = e^p* for the 
brane at y = 1/2, Eqs. @ and (|8|) read 
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Had we introduced a quadratic term in y in the ansatz for the BD field in Eq. (24), we would end up with a third 
constraint 

[0']o = -[0']l/2. (31) 

We must note that this constraint is formally equivalent to ( P7| ) and ( p8|) in the sense that the three constraints require 
the jump in the first derivative of a given quantity across one brane to have the same absolute value but an opposite 
sign as the same jump across the second brane. Since we have only two quantities to fix, namely p* and p*, this 
extra relation would possibly be verified only for a discrete set of values of lo and/or e. Therefore by discarding the 
quadratic term in we are automatically free from this problem. In any case our ansatz can always be interpreted 
as a power series expansion in y of a more complex solution in which case it is valid for y small enough. Since we are 
mainly interested in the cosmology of the brane at y = we are well within the domain where it is safe to consider 
only the lower order terms in the power series. 

Although we are not going to discuss further the problem of constraints in this paper, we suggest that a much 
more elegant solution could possibly be found with a varying uj in the bulk. This could in principle be implemented 
in the context of the so called scalar-tensor theories where uj = ^^-"(0). The third constraint would then fix the values 
of w at y = 1/2 as a function of uj on the brane at y = 0, therefore partially constraining the functional form of 
Lu = uj{(j)). Another possibility which would make the constraints unnecessary is a model with a single brane. In that 
case, Eq. (^) would hold without any further ado. We must point out however that the presence of a second brane is 
favored for several reasons (see Refs. 0,0). As a final remark, let us note that constraints of this type are expected 
to appear in models with compact extra dimensions |q| . 



V. EXACT COSMOLOGICAL SOLUTIONS 

We will now obtain some simple solutions for a and which will allow us to discuss cosmology in our setup. For 
what follows we will discard the part of our ansatz in Eq. (g^) (and similarly for n) which is proportional to y^. In 
this context this means that we will no longer have to worry about the constraints discussed in the previous section. 
As before we assume the extra space like dimension is stabilized (6 is constant) and choose time such that ng = 1. 

From Eq. (|26| ) and using the jump conditions Eqs. (|l7|-|l9|) we obtain the Friedmann equation which reads 
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Using our ansatz in the BD field equation, we obtain again Eq. (E2|), as we should expect. Together with the 



conservation equation (gO|), Eqs. (g2|) and ( |32[ ) form a complete system which will allow us to study the cosmology 
on the brane. A careful look at the above mentioned equations suggests that once the equation of state on the brane 
has been specified, it is possible to obtain a power law solution to the scale factor, BD field and the matter content 
on the brane as a function of time. 

Here we consider only the particular cases of radiation (p — p/i) and dust (p = 0) as well as the vacuum energy 



case (p = — p). For a radiation dominated Universe we obtain 0o — const, from Eq. (22), which also holds for the 
usual BD solution in four dimensions. It is noticeable that in this case the solution is exactly the same as in the case 
of general relativity (see Ref. |g]): 

flo oc T""^' ; Pa (x a^ ; ipo — const . (33) 

For the matter dominated era the solution can be obtained from Eqs. (pO[), ( [2^ ) and (p2|), and is given by 

3^ + 3+%/3v'2^ + 3 l+y3\/2^ + 3 

aa (X T 5;:;tt2 ; p^ oc a^-^ ; 0o oc r ^=^ , (34) 

where we have taken the growing mode for the scale factor and we have restricted to large tu by truncating terms 
0{l/u;^). It is important to notice that for large w the solutions approach the brane general relativity case a oc r^/^. 
If the solution for the scale factor in the 4 dimensional BD case is written as a ex r" then in our model the solution 
takes the form a oc r"'^+^^' ^'. In our situation however the expansion during the dust era is faster compared with 
that in the general relativity case. This will have some implications during the rapid oscillations of the scalar field in 
our brane after inflation and during preheating when the expansion mimics that of a pressureless fluid. 

It is important to mention here that our results are valid only prior to nucleosynthesis and at late times the 
Friedmann equation should converge to the standard case with H oc ^ in order to satisfy the stringent constraints 
coming from big bang nucleosynthesis. For this to happen one may need to take the brane tension into account [pj. 
Denoting the brane tension by /5a, the total energy density in the brane will be given by ptot — PA + P and the 
square of this quantity will include the term 2p\p linear in p. For pA > IMeV one can ensure compatibility with 
nucleosynthesis constraints. In this way for energy scales larger than the brane tension we have the non-standard 
behaviour, while as the universe expands and cools down we recover the standard behaviour. 

Finally, we discuss the cosmology of the vacuum energy dominated stage. From the energy conservation equation 
we see that false vacuum inflation {p =constant) is possible in this case for large enough values of w. The scale factor 
and BD field are given by 

uq^t) oc r^e (f)Q[T) oc T . (35) 

The rate of expansion during inflation in the brane is therefore slower when compared with the 4 dimensional BD 
case a oc t^+'^/^. As for the BD field, in the 4 dimensional case oc r^ and the evolution is also slower in the brane 
model. 

Note that if the energy conservation equation (po| ) was not satisfied as in the models studied in P2|-p4| then false 
vacuum inflation would not be possible because due to matter leaking from the brane, even for the case p = —p it 
would not possible to keep the energy density constant anymore. 

An interesting possibility which we have not analyzed here is the case where the dynamics in the brane is driven 
by the evolution of the fifth dimension. A more detailed study of the cosmology in the model described in this paper, 
including an analysis of the primordial perturbation spectra produced by infiation in our setup will be discussed 
elsewhere. 



VI. CONCLUSIONS 

We have derived the effective four dimensional field equations for the brane cosmology in presence of the BD field 
in five dimensions. In our scenario, the conservation equation for the matter field stuck to the brane still holds in 
spite of the BD field present in the bulk. This is due to the fact that our calculation was done in the Jordan frame. 
Some global solutions with two branes such that the fifth dimension has an orbifold symmetry were discussed. Finally 
we obtained the modified Friedmann equation in the brane and cosmological solutions in the BD brane setup were 
discussed. 
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